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ABSTRACT 

We present results from JCMT sub-mm observations of sources selected from the ISO 
FIRBACK (Far-IR BACKground) survey, along with UKIRT near-IR imaging of a 
sub-sample. This gives valuable insight into the brightest ~10% of galaxies which 
contribute to the Cosmic Infrared Background (CIB). We estimate the photometric 
redshifts and luminosities of these sources by fitting their Spectral Energy Distribu- 
tions (SEDs). The data appear to show a bimodal galaxy distribution, with normal 
star-forming galaxies at z ~ 0, and a much more luminous population at z ~ 0.4-0.9. 
These are similar to the ultraluminous infrared galaxies which are found to evolve 
rapidly with redshift in other surveys. The detectability threshold of FIRBACK bi- 
ases the sample away from much higher redshift (z > 1.5) objects. Nevertheless, the 
handful of z ~ 0.5 sources which we identify are likely to be the low- 2 counterparts 
of the typically higher- z sources found in blank field sub-mm observations. This sub- 
sample, being much more nearby than the average SCUBA galaxies, has the virtue of 
being relatively easy to study in the optical. Hence their detailed investigation could 
help elucidate the nature of the sub-mm bright galaxies. 

Key words: infrared: galaxies - submillimetre. 



1 INTRODUCTION 

What makes up the Cosmic Infrared Background (CIB) 
detected from the COHE-FIRAS data (Puget et al. 1996; 
Fixsen et al. 1998; Lagache et al. 1999; Lagache et al. 2000; 
Hauser et al. 1998; Dwek et al. 1998; Finkbeiner et al. 
2000)? This remains an open question, and details of galaxy 
types, their redshift distribution, and how they appear in 
other wavebands, remain sketchy. The FIRBACK (Far-IR 
BACKground) survey (Puget et al. 1999; Dole et al. 2001) 
addressed this question by performing some of the deep- 
est blank field ISO surveys at 170 (jm, near the peak of 
that radiation. About 200 sources were detected above 3 a 
(=135 mjy) accounting for about 7% of the nominal back- 
ground value. 

In general far-IR sources such as the FIRBACK ones 
sample the low-to-moderate redshift regime, and thus pro- 
vide a link between the local Universe and high-z sources, 
such as the SCUBA-bright 'blank-sky' population (see Blain 
et al. 2002 for a review). Understanding the nature of these 
sources, their emission mechanisms, and their dust proper- 



ties is crucial to our understanding of galaxy formation and 
evolution from high redshift until today. This in turn informs 
us about the cosmic background, as well as the nuclear ac- 
tivities, star formation distributions, and the role of dust 
obscuration in star formation through a large fraction of the 
history of the Universe. 

In order to better understand the sources detected by FIR- 
BACK we have been carrying out follow-up observations 
with SCUBA at 450 (im and 850 (im. Detection by ISO at 
170 /im (5'i7o>135 mjy), means that a strong bias away from 
high-z objects is present, although we still expect to detect 
objects out to z ~ 1. The FIRBACK galaxies represent the 
brightest contributors to the CIB, and are a different se- 
lection than typical SCUBA galaxies. 'Blank-sky' sub-mm 
bright galaxies, although so far accounting for up to 50% of 
the sub-mm background (e.g. Cowie, Barger, & Kneib 2002), 
make up an insignificant fraction of the total CIB. The com- 
bination of far-IR and sub-mm observations is thus very 
powerful in establishing a link between high-z dusty star- 
bursts and their local counterparts. 

Observationally, without redshifts, it has been difficult 
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to distinguish between cooler local starbursts and warmer, 
more luminous sources at higher z. This is because the spec- 
tral energy distributions (SEDs), for a fixed emissivity index 
/3, are degenerate in the parameter combination (1 + z)/Td. 
There are additional complications of course caused by vari- 
ations in j3, Td, and luminosity, changing the shape of the 
SED. Fundamentally, our understanding of dust in extra- 
galactic sources, its properties, and interaction with the ra- 
diation field is poor, which is a major impediment in our 
interpretation of the observational evidence. 

The only way to approach these issues is by detailed 
multiwavelength studies of samples representing key ele- 
ments in the above puzzle. 

The data we present in this paper constitute a far-IR 
selected sample. This sample is now large enough, and with 
wide enough wavelength coverage, including near-IR, far-IR, 
sub- mm and radio observations, to be able to tackle some of 
these issues. We do this through a combination of direct SED 
fitting, statistical analysis, consistency with other observa- 
tions and comparison with model predictions, trying to use 
the minumum number of a priori assumptions. In essence, 
we start with the assumption that our sources represent a 
uniform population, which we demonstrate are most likely 
local starforming galaxies. This model breaks down for a 
handful of our sources, which is probably due to their being 
at somewhat higher redshifts (and more luminous) than the 
rest of the sample. Since our sample is representative of FIR- 
BACK as a whole, we have thus improved our knowledge of 
the nature of the brightest sources making up the CIB. 

We focus largely on the conclusions deriving from the 
sub-mm observations. Discussion of the individual proper- 
ties of FIRBACK galaxies in general, including issues re- 
lating to the identification of the sources will be discussed 
elsewhere (Lagache et al. 2003). 

Throughout the paper we assume a flat Universe with 
#o=75 km s _1 Mpc -1 , n M =0.3, and fi A =0.7. 



2 THE DATA 
2.1 Sample 

Targets were selected from the FIRBACK (> 3a) catalogue 
(Dole et al. 2001) in the ELAIS Nl field. Confusion remains a 
major issue with a beam size ~ 90 arcsec. In our selection we 
aimed at a uniform sample that would be as representative 
of the FIRBACK population as possible. Our main selection 
criterion has been the availability of radio detections (Ciliegi 
et al. 1999) inside the beam (required for JCMT pointing). 

We exclude FIRBACK sources with radio-bright coun- 
terparts which are rare and not representative of the typical 
FIRBACK source, as well as sources with several radio de- 
tections per beam *. This left us with 41 possible sources 
to draw from - we have followed up 30 (and additionally 
one N2 source) - listed in Table 1. Due to this broadness 
of selection criteria we believe our sample represents a fair 
cross-section of the FIRBACK population, rather than fo- 
cusing on a specific sub-population. This can be seen clearly 



* This criterion was not applied too strictly and a few sources 
(Nl-008, Nl-029, Nl-041) do in fact have two radio sources inside 
the beam. 
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Figure 1. Distribution of the S170/1S1.4GHZ ratio. The solid his- 
togram represents the 41 sources from the FIRBACK catalogue 
which meet our selection criteria. The dashed histogram is our 
sub-sample of 31 targets. The dotted line shows our high-z can- 
didates (discussed in Section 4). 



in Fig. 1, where there is close agreement between the distri- 
bution of S'i7o/5'i.4GHz for all 41 possible sources with that 
of the 30 sources in our sample. 



2.2 SCUBA sub-mm observations 

The observations presented here were taken with the Sub- 
mm Common User Bolometer Array - SCUBA (Holland et 
al. 1999) instrument on the James Clerk Maxwell Telescope 
(JCMT) in March 1999 and in March and May 2001. In or- 
der to avoid biasing our data, we attempted to observe each 
source until a predetermined rms (~ 1.5 mjy) was reached, 
irrespective of whether the source appeared to be a possi- 
ble detection or not. The March 2001 data were taken in 
exceptional grade 1 weather (T225 ~ 0.04, and as low as 
0.02), whereas the 1999 data were taken in merely 'good' 
weather (T225 ~ 0.07). Throughout our observations we used 
the 2-bolometer chopping mode. This involves chopping in 
array coordinates in order to always align one negative beam 
exactly with a specific off-centre bolometer with the other 
being partially aligned. Thus the negative beams ca n be 
folded in, improving the rms by a factor of up to \/2/3. 
Unfortunately, for our 2001 run, software problems with the 
new telescope control system resulted in SCUBA not prop- 
erly chopping onto another bolometer, making the negative 
beams unrecoverable. Thus the 2001 data presented here are 
from the central bolometer only, whereas the 1999 data have 
the negative beams folded in. 

The data were reduced using the SURF package (Jenness et 
al. 1998), and also with custom- written code. The extinction 
correction was performed using sky-dip observations when- 
ever available, and with a derived optical depth from the 
Tcso — Tscuba relations (Archibald et al. 2000) otherwise. The 
sky was subtracted as a mean level for the entire array at 
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each jiggle, excluding the signal bolometers as well as any 
bad bolometers, as revealed by excessive noise levels. A cal- 
ibration uncertainty of a maximum of 20% exists (based on 
differences among the nightly calibration values), however 
it has little effect on the signal-to-noise ratios, especially for 
the 2001 data where each source was typically observed in 
a single night. The pointing uncertainty was ~ 2 arcsec. The 
SCUBA 850 /xm beam FWHM is ~ 15 arcsec, whereas the 
450 fim beam is ~ 8 arcsec. 

The derived 450 fim and 850 fim flux densities and associ- 
ated 1 a errors are listed in Table 1. 



2.3 UKIRT near-IR observations 

Most of our SCUBA detected F1RBACK sources could not 
be identified in DPOSS images. As the sources are ex- 
pected to be extinguished by dust and therefore have red 
spectra, we obtained deep observations in the K-band at 
UKIRT, using the Fast Track Imager (UFTI) for maximum 
sensitivity to obscured components. The small UFTI field 
(50 arcsecx 50 arcsec) was centred on the source positions 
taken from the radio/SCUBA identifications. Each source 
was imaged for a total of 1800 s, with individual exposures 
of 60 s each, reaching a limiting magnitude in a 2 arcsec di- 
ameter aperture of K = 20.4 (5a). The fast tip/tilt, adap- 
tively corrected imaging resulted in seeing better than me- 
dian conditions, at ~0.4 arcsec FWHM. Data were reduced 
using the Starlink, UKIRT/UFTI image processing tools un- 
der the ORACDR environment (Bridger et al. 2000). We wrote 
custom ORACDR scripts to optimize point source sensitivity 
in our essentially blank field observations, creating flat fields 
from each 9-point dither, and high signal-to-noise thermal 
background images from 60 minutes of data centered around 
the observing period of each target. Fig. 2 shows the avail- 
able UKIRT images. For the rest of our sample we use data 
from the 2MASS catalogue (for sources with K ~ 14-15 we 
estimated their magnitudes directly from the calibrated cat- 
alogue images via aperture photometry in GAIA). 



2.4 Far-IR and radio fluxes 

IRAS 100 (jm fluxes were obtained using the xscanpi facil- 
ity. We quote them here for the sake of completeness, as well 
as to help in comparison with local IRAS galaxies (which 
are typically >Uy). The errors quoted are purely statisti- 
cal; however, at this faint level systematic errors (including, 
among other things IR cirrus and mapping artefacts) dom- 
inate the flux estimates, to the point of making them of 
little statistical use. ISOPHOT 90 fim detections also exist 
for a few of the objects, as well as mid-IR ISOCAM data. 
We do not use these data as they are too sparse (typically 
<5 sources) and provide little constraint on the thermal far- 
IR/sub-mm, which is our main focus here. 
The radio fluxes at 1.4 GHz given in Table 1 are from a VLA 
survey of the field (Ciliegi et al. 1999), and errors are 1 a. 



3 RESULTS & ANALYSIS 

3.1 Assembling the multiwavelength data 

The 1999 data presented here were previously discussed in 
Scott et al. (2000)*. However we re-reduced the old data 
concurrently with the 2001 data in order to ensure unifor- 
mity, especially as an upgraded version of SURF and a new 
custom-written code were used. We confirm all previously 
reported detections (i.e. targets having SNR>3). The 2001 
data have three unambiguous detections at 450 /im (N 1-004, 
N 1-024, N 1-078). This high detection rate at a difficult band 
is due to the exceptional atmospheric conditions during our 
observing run, as well as the far superior performance of the 
new wide-band filter. In addition there are three new detec- 
tions at 850 /im (Nl-001, Nl-059, N-078). Note that the few 
arcsecond pointing uncertainty (see Section 2.2) has only a 
small effect on the long wavelength data, but may be a signif- 
icant reason for the apparently missing 450 fim flux (where 
the beam FWHM is only ~ 8 arcsec) in sources where one 
would expect to find some (e.g. Nl-059). Our results for the 
entire sample are presented in Table 1, which, in addition to 
the sub-mm data, includes near-IR, far-IR, and radio data. 
The sample as a whole is strongly detected, with aver- 
age fluxes of (5 , 85o)=2.6±0.2 mjy and (S , 4 5o)=16.7±1.9 mjy. 
This is indicative of the fact that although most sources do 
not clear the 3a detection threshold, there is useful statisti- 
cal information in the flux data (far in excess of what would 
be expected for mere sky fluctuations). Thus we treat these 
fluxes directly, along with their uncertainties, rather than 
just using the upper limits. Here we focus on the general 
properties of the sample which are revealed by the addition 
of the sub-mm and near-IR data. Physical properties for in- 
dividual sources will be explored elsewhere (Lagache et al. 
2003). 

Two possible sources of uncertainty in comparing the far- 
IR and SCUBA fluxes are: 1) the different sizes of the 
ISO(about 90 arcsec) and SCUBA (about 15 arcsec) beams, 
and thus the possibility of multiple sources lying within the 
beam - this is also complicated by the issue of clustering, 
and possibly lensing; and 2) the relative proximity of the 
bulk of our sources and thus the possibility of at least a few 
being extended beyond the SCUBA beam (see for exam- 
ple Nl-004 and Nl-007 in Fig. 2). This is especially severe 
for the ~ 8 arcsec 450 fim beam. In general, however, images 
at other wavelengths, are insufficient to determine precisely 
how much flux is being lost in this way, as it is not clear 
how concentrated the sub-mm emission would be, or how it 
is distributed in comparison with other wavebands. In any 
case, the sub-mm flux estimates are likely to be less accu- 
rate in photometry mode if the source size is on the order 
the beam size. These effects, of course, become more impor- 
tant with increasing Si 70. 

Although the low sub-mm flux of most of our sources is 
consistent with local star-forming galaxies, there is the addi- 
tional concern about the possibility of cirrus contamination 
of the FIRBACK sample. Previous studies have discussed 



* Note that we use the naming scheme of Dole et al. (2001), 
which differs from the earlier convention used. In particular Nl- 
038, Nl-061, and Nl-063 from Scott et al. (2000) correspond to 
the new Nl-040, Nl-048, and Nl-064. 
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Figure 2. The available UKIRT iiT-band images. The images are centered on the radio position, and are about 15xl5arcsec or roughly 
a SCUBA 850 (im beam size. Exceptions are Nl-004 and Nl-007, where the white circles have a diameter of ~15 arcsec. 



this issue extensively (Lagache & Puget 2000; Puget et al. 
1999), concluding that the FIRBACK sources are unlikely 
to suffer from significant cirrus contamination. The Nl field 
was chosen for this survey due to its low cirrus emission. We 
can estimate, using the standard formula (Helou & Beich- 
man 1990), that the cirrus variance is ~10 times less than 
the detection cutoff of the FIRBACK survey. Due to the 
non-Gaussianity of the cirrus fluctuations (Gautier et al. 
1992) some outliers are still possible, but these would have 
to reside in very cold cores which have not been revealed in 
deep searches for CO emission. Moreover, the far-IR slope 
of cirrus knots is known to be much steeper than that of 
sources. A study, whose source selection is similar to FIR- 
BACK's (Juvela et al. 2000), makes use of this in concluding 
that cirrus contamination is highly unlikely. In agreement 
with that study, our sources show a fairly flat Sioo / £170 ra- 
tio. Finally, all our FIRBACK sources have a radio counter- 
part, further reducing the chance that cirrus contamination 
is relevant. 



3.2 Linear Correlations 

Fig. 3 shows three illustrative projections for our data, using 
the 170 /im, 450 /im, 850 /im and radio fluxes. For each plot a 
fit is made to all points using a straight line passing through 
the origin, i.e. y — mx. Then the > 2a outliers are excluded, 
a new fit is made and is plotted along with the rms scatter. 
Table 2 shows the results. Since in most cases, the errors 
in the two directions are of comparable size, rather than 
use the conventional 1-D x 2 , we minimize a statistic which 
combines errors in both directions: 



2 (y l -mx i ) 2 

2 — 1 y 

where a x and a y are the two errors. In general, it appears 
that the typical uncertainties of our sub-mm data points are 
larger than the 1 a scatter of the fits. We test this explicitly 
for the S450 vs. Sgso plot via Monte Carlo simulations of the 
data, thereby estimating the probability of the X2D obtained 
to be ~ 2%. This implies that our errors are overestimated 
by a factor of ~ \/2 or else there is some as yet unidentified 
source of correlation in the errors. We will, however, be con- 
servative and leave the errors as they are, since we cannot 
properly account for the source of this discrepancy. 
After this cautionary aside, we return to Fig. 3. It is orga- 
nized such that the order roughly follows the main features 
of the long wavelength SED: the top plot largely tracks the 
location of the thermal peak; the middle corresponds to the 
sub-mm slope; and the bottom plot traces the trough be- 
tween the thermal and non-thermal emission. 
First, we will concentrate on the Ssso vs. Si 70 plot (Fig. 3a). 
We immediately notice that about 5 sources (the circled 
crosses) occupy a locus > 2a away from the best-fit line. 
Assuming a grey-body model, these sources are either at 
higher redshift or lower temperature than the rest of the 
FIRBACK 170 /im sources (or a combination of both). 
To understand which, we turn to the 5*450 vs. Ssso plot 
(Fig. 3b). Here we are mainly exploring the slope of the 
spectrum in the sub-mm. What we notice in this plot is that 
there is not a population of outliers as there was in the 5*850 
vs. S170 plot (the only outliers are Nl-001, Nl-034, and Nl- 
059 which all have sub-mm slopes which alone cannot be fit 
by any sensible dust/redshift combination and thus we as- 
sume suffer from some systematic effect such as discussed in 
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Table 1. Multiwavclength data for our sample. All errors are 1<t estimates. 



Source 




Sioo 6 




5 450 d 


5 8 50 d 


■5l.4GHz e 






mag 


mjy 


mJy 


mJy 


mJy 


mJy 


M- 


-001 


12.4±0.1 


430i87 


597i72 


-3.0il4.0 


6.1il.6 


0.74i0.23 


Nl- 


-002 


12.7i0.1 


340il21 


544i69 


14.4il2.4 


4.4il.l 


0.64i0.04 


Nl- 


-004 


12.4i0.0t 


300i73 


391i58 


32.5i7.2 


3.6il.4 


0.88i0.13 


Nl- 


-007 


13.2i0.lt 


480i73 


338i54 


23.4i8.1 


4.4il.6 


1.04i0.12 


M- 


-008 


14.2i0.lt 


160i73 


335i54 


25.9il4.0 


1.9il.l 


2.98i0.04 


Nl- 


-009 


12.0i0.0 


310i58 


313i52 


10.6i7.6 


3.5il.5 


1.15i0.11 


Nl- 


-010 


13.0i0.0 


360i99 


309i52 


15.2il0.8 


1.8il.4 


1.05i0.20 


Nl- 


-012 


13.9i0.2 


320il22 


302i51 


9.2il0.0 


1.5il.6 


0.31i0.07 


Nl- 


-013 


16.8i0.lt 


350i75 


294i51 


18.8i9.9 


0.0il.5 


0.52i0.15 


Nl- 


-015 


14.8i0.1 


230i41 


294i51 


-3.4±7.7 


1.4il.6 


0.52i0.07 


Nl- 


-016 


13.2i0.1 


360i92 


289i50 


34.8il6.7 


1.5il.2 


1.55i0.13 


Nl- 


■024 


14.2i0.0t 


<147 


266i49 


32.3i7.5 


2.9il.3 


0.75i0.02 


Nl- 


-029 


14.3i0.1 


340i83 


229i46 


20.0il4.2 


0.5il.7 


0.69i0.05 


Nl- 


■031 


13.5i0.1 


110±78 


225i46 


9.2il3.2 


1.9il.l 


0.43i0.06 


Nl- 


-032 


18.5-19. 5t 


220i53 


224i46 


14.9i7.7 


1.3il.4 


0.21i0.05 


Nl- 


■034 


19.3i0.8t 


270i71 


221i46 


95.1i27.5 


1.3il.3 


0.33i0.07 


Nl- 


■039 


15.8i0.2 


230i 


205i44 


10.9i86.8 


-0.1i2.3 


0.58i0.07 


Nl- 


■040 


19.4i0.6t 


< 153 


205i44 


29.2i20.5 


5.4il.l 


0.33i0.03 


Nl- 


-041 


14.7i0.1 


150i55 


204i44 


20.4il56.7 


-0.1i2.5 


0.76i0.06 


Nl- 


■045 


14.4i0.0t 


< 165 


198i44 


15.3i8.3 


3.0il.4 


0.43i0.06 


Nl- 


-048 


19.3i0.9t 


< 165 


192i44 


15.5ill.6 


4.2il.l 


0.37i0.05 


Nl- 


■056 


15.9i0.2 


160i71 


179i43 


-8.6i8.4 


0.0il.6 


0.24i0.02 


Nl- 


■059 


20.4i0.9t 


230i66 


175i42 


6.9i34.5 


6.4il.9 


0.60i0.06 


Nl- 


■064 


18.2i0.3t 


260i79 


166i42 


35.2il3.9 


5.1il.2 


0.23i0.04 


Nl- 


■068 


15.3i0.lt 


320i95 


165i42 


15.1±7.6 


2.2il.4 


0.44i0.05 


Nl- 


-077 


15.5i0.2 


200i89 


159i41 


5.9i7.3 


l.lil.3 


0.40i0.10 


Nl- 


-078 


18.0i0.4t 


240i63 


158i41 


35.2i8.7 


5.7il.3 


0.24i0.04 


Nl- 


-083 


15.4i0.2 


< 195 


150i41 


16.2il6.0 


0.7il.2 


0.55i0.03 


Nl- 


■101 


15.2i0.2 


210i73 


136i40 


19.8i7.5 


0.9il.5 


0.39i0.05 


Nl- 


■153 


15.5i0.2 


140i58 


103i37 


9.6il5.3 


-0.2il.O 


0.24i0.03 


N2- 


■013 




310i75 


244i53 


23.5il5.9 


3.5il.4 


0.30i0.07 



a These were obtained from our UKIRT sample (marked by t, see also Fig. 2) when 

available and from the 2MASS catalogue otherwise. Those fainter than K ~ 14 were 

estimated directly from the 2MASS images. 

6 IRAS 100 /im fluxes obtained using XSCANPI (see text). 

c ISOPHOT 170 ^m data from Dole ct al. (2001). 

d SCUBA 450 ^m and 850 (im fluxes from this work. 

e VLA 21cm fluxes from Ciliegi et al. (1999). 



section 3.1). This common distribution implies that a single 
[f3,T/(l + 2)] combination describes the sample reasonably 
well. The lack of the outliers, from the Ssso vs - <5i70 plot 
(Fig. 3a), means that the 5 outliers' (Nl-040, Nl-048, Nl- 
064, Nl-059, and Nl-078) location in Fig. 3a is most likely 
due to having somewhat higher 2 rather than a significantly 
different SED shape. In addition, their redshifts can only 
be as high as 2 ~1, corresponding roughly to the shift be- 
tween 170 /im and 450 /im (a higher redshift, much warmer 
population would have to be very finely tuned to still fall 
on the same distribution, which seems improbable) . This in 
turn roughly constrains the temperature of the sources (by 
Wien's law) to less than about 40 K. 

We now turn to Fig. 3c, the 850 pm/1.4 GHz vs. 850 /tm 
comparison, where the previous outliers return again (ex- 
cept for Nl-059 which has unusually high radio emission, 
discussed later), confirming that indeed redshift, rather than 
dust properties is the main difference between them and the 
rest of the sample. 



3.3 Sub-mm/Radio redshifts 

As a tracer of the trough between the thermal and non- 
thermal emission of galaxies, Sssq/Siaghz is commonly used 
as a redshift indicator (Carilli & Yun 2000). It is, however, 
degenerate in dust properties for galaxies cooler than about 
60 K (Blain et al. 2002) . Since we are clearly in that regime 
(see previous section), we have investigated using this re- 
lation calibrated on samples with different selections. The 
first parameterization we used was the Carilli & Yun re- 
lation (Carilli & Yun 2000, CY hereafter) which is based 
on JHAS-selected and somewhat more radio-loud galaxies, 
with likely a higher fraction of AGN. The second parame- 
terization we used was that from Dunne, Clements & Eales 
(Dunne, Clements & Eales 2000, DCE hereafter) which is 
based on the SLUGS (SCUBA Local Universe Galaxy Sur- 
vey) sample (Dunne & Eales 2001) which is essentially an 
ffiAS-bright local selection (unlike ours). 
In general, such relations provide too weak a constraint on 
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Figure 3. Multiwavclcngth projections of our data. The solid lines 
are the y = mx fits obtained by minimizing x|d equation (1). The 
dashed lines are the ±1<t scatter of the points around the lines, 
with the fit parameters given in Table 2. The error bars shown 
are representative for our data. The 5 outliers in panel (a) are 
indicated with the same symbols in the other two panels. 



Table 2. Results for the linear fits to the data 

Relation m rms \ 2 

850 fim vs. 170 (im 0.01 1.23 18.07 

450 fim vs. 850 7.05 11.74 13.95 

850 Mm/1.4 GHz vs. 850 Atm 1.82 2.46 12.19 



individual redshifts locally, due to the large scatter in intrin- 
sic galaxy properties. However, they still confirm our prior 
selection of Nl-040, Nl-048, Nl-064, Nl-059, and Nl-078 as 
likely being at somewhat higher redshifts than the rest of 
the galaxies in our sample (based on both relations giving a 
redshift > 0.4 for all of these sources). The only other source 
which satisfies this is N2-013. We return to it in section 3.6, 
when we also discuss different redshift estimators. 



3.4 Sub-mm vs. near-IR 

Here we examine the correlation between the K magnitudes 
and 850 fim fluxes of our sources. The 5*850 flux density by 
itself is not a good redshift indicator due to its fc-correction 
behaviour. It is however a good luminosity tracer. The K 
magnitude (in the rest frame) is also a luminosity indicator, 
as it is 10 times less obscured than the optical. The Ssso- 
to-K magnitude relation can be used as a redshift indicator 
(Barger et al. 1999; Dannerbauer et al. 2002) because the 
rest-frame shorter wavelengths (that are much more dust 
obscured) move with increasing redshift into the observer- 
frame near-IR. Thus the K magnitude for a given sub-mm 
flux is dependent on both redshift and dust obscuration. 
Fig. 4 plots the flux at 850 fim against K magnitude, where 
the higher- z sources clearly populate a different locus from 
the nearby galaxies. They are roughly 3 a removed from the 
best-fit relation for the other sources, with a distinct gap be- 
tween the two groups (apart from a couple of sources which 
we discuss later). The gap is more pronounced than in Fig. 3, 
as here the (1 + z)/T d degeneracy is somewhat broken. 

A redshift relation based solely on the K magnitude is 
bound to be degenerate in some other galaxy properties 
(amount of dust, luminosity etc.), and can only work for 
a very homogeneous sample of galaxies (see e.g. Willot et 
al. 2002). The addition of sub-mm flux improves this rela- 
tion, since the if-band is dust absorption attenuated (specif- 
ically for distant, luminous galaxies), whereas the sub-mm 
flux arises from dust emission, and thus the combination of 
the two will somewhat break the dust degeneracy. This diag- 
nostic of both the absorption and emission spectrum allows 
for a more robust redshift indicator over different galaxy 
types. What is robustly clear is that, in general, objects 
which are detected at 850 pm and are faint at AT-band are 
at higher redshift. In Fig. 4 we find the same handful of 
outliers as in Sections 3.2 and 3.3. 



3.5 SED fits 

We now fit single grey-body SEDs to the 170 /im, 450 /im, 
and 850 ^m fluxes of each source (see Fig. 5). We assume 
optically-thin sources; the effect of including a non- negligible 
r in the fits is to suppress the peak with respect to the 
Rayleigh- Jeans tail, such that the best-fit dust temperature 
inferred will be 10-20% higher than otherwise (Blain et al. 
2002). This is therefore not an important effect for sub-mm 
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spectra dominated by single temperature dust emission (to 
the precision of our fits of three points). In order to avoid 
stretching these assumptions too much, we only use the 
170 fim, 450 /im, and 850 /tm fluxes in the fits, where the 
SED is dominated by the coldest significant dust compo- 
nent. 

We fit for two parameters - the overall normalization, which 
gives the luminosity (if the redshift is known), and the wave- 
length shift, which is proportional to (1 + z)/Ta. The best 
fit x 2 arm fit parameters were obtained with the pikaia ge- 
netic algorithm (Charbonneau 1995). With only 3 points 
for each galaxy, we could not also fit for the emissivity in- 
dex f3, and therefore it was held constant. However, we in- 
vestigated how different values of /3 (1.3, 1.5, 1.7, and 1.9) 
would affect the quality of fits for the sample as a whole. 
Values of P = 1.5 — 1.7 are far better than either smaller or 
larger values, for which a third of the sample has a x 2 > 2. 
The sources, with a poor fit in all cases are Nl-001, Nl-002, 
Nl-015, Nl-034, and Nl-056. For the lowest /3 tested, 5 ad- 
ditional sources show a poor fit: Nl-008, Nl-009, Nl-012, 
Nl-013, Nl-016 (note that these are all among the brighter 
ISO sources in our sample) . On the other hand, the higher, 
P — 1.9, value, provides a poor fit for all the sources we 
singled out in earlier sub-sections as potentially being at 
higher redshifts (see Fig. 5). For local ULIGs, previous stud- 
ies (Dunne & Eales 2001; Klaas et al. 2001) already showed 
this trend, with lower /3 providing a better fit when single 
temperature grey-bodies are used, whereas generally a f3 ~ 2 
is inferred for multi-temperature fits. 

Notice that in Fig. 5 we also show the 100 /tm point, al- 
though it is not included in the fits. Fitting to all four points 
results in poor fits either on the shorter or longer wavelength 
ends, depending on the choice of /3. To fit this wider range 
presumably requires a multi-component model, but the data 
do not warrant a detailed investigation of more complicated 
SEDs. 

It thus appears that it is reasonable to model our data as a 
single temperature grey body in the limited spectral range 
considered (i.e. 170 ^im^850 /im), provided we explore the 
effects of varying the parameters within the accepted range 
we found here. Of course this model will be incorrect in detail 
for the entire thermal spectrum, as has been shown previ- 
ously by several authors (Dunne & Eales 2001; Blain et al. 
2002; Stickel et al. 2000; Klaas et al. 2001). In general, any 
grey-body model, even with multiple components, remains 
only an approximation to the underlying far more complex 
dust properties. We also caution that the dust properties 
are merely phenomenological, corresponding to emission- 
weighted averages, and depending on the wavelength range 
chosen. Care should therefore be taken in intcrpretting any 
dust temperatures physically. 

3.6 Redshifts, and Luminosities 

Now we wish to examine some of the physical characteris- 
tics of the non-zero redshift sources in our sample (because 
they are the best constrained through their sub-mm fluxes) . 
Since, fundamentally, the dust spectrum/distance degener- 
acy is still present, derivations of properties such as luminos- 
ity have an intrinsic uncertainty. However, we can provide 
physical parameter estimates for what we consider to be 
reasonable ranges of redshifts (for the likely higher redshift 




K magnitude 



Figure 4. SgSO vs K magnitude for our sample (crosses, with the 
5 outliers also encircled). The filled square and ii"-band upper lim- 
its are from the SCUBA lens survey (Smail et al. 2001), while the 
pentagons are from the UK 8 mjy survey (Ivison et al. 2002). Our 
high-z candidates populate a similar region to these other SCUBA 
survey sources. The solid line is a log-log fit to the crosses only, 
with the dashed lines being the ±1<T scatter. The error bar in the 
lower right-hand corner is a representative one for our measure- 
ments. 

sources only) in order to obtain a handle on the physical 
nature of our sources for comparison purposes. 
A few of our optically-fainter sources have available spec- 
troscopic redshifts obtained from Keck and Palomar spec- 
troscopy. These are Nl-045(z=0.25); Nl-068(«=0.22); Nl- 
039(2=0.27); Nl-008(z=0.27); Nl-040(,z=0.45); and Nl- 
064(2=0.91). The details of the spectroscopic observations, 
and the obtained redshifts are discussed elsewhere (Chap- 
man et al. 2003; Chapman et al. 2002). These provide valu- 
able checks of our redshift estimates. However, the spec- 
troscopic redshifts only exist for a fraction of our sample, 
having strong selection effects. Therefore, we do not make 
more direct use of them as we wish to treat our sample in a 
uniform way. The details of the spectroscopic observations, 
and the obtained redshifts are discussed elsewhere (Chap- 
man et al. 2003). 

For definitiveness we focus on the (/3, T d ) combinations 
which agree with the SED fit results (see section 3.5), with 
the redshifts estimated from the sub-mm/radio relations 
(section 3.3), and with expectations from other studies (see 
section 4). These are (1.7, 30 K), and (1.5, 40 K). Table 3 
shows the redshift ranges, along with luminosities^, and cor- 
respoinding distances for the higher redshift candidates in 
our sample. The ranges illustrate how the values vary with 
input parameters, and also give an indication of the spread 
of reasonable values. The sources in the upper half of the 

T We solve for the total far-IR luminosity simply by integrating 
our fitted grey-body SEDs. The key mathematical step is the 
integral x s ~ 1 (e x - l)" 1 ^ = T(s)C(s). 
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Table 3. Fit results for all the higher- z sources" 
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Figure 5. The SED fits for a fraction of our sample where the fits 
suggest a non-zero redshift. The axes scales are linear with the 
x, {/-ranges in all panels being the same (shown in the labels), and 
the sources are arranged by decreasing 170 fim flux. The IRAS 
100 jim points included in the plots are not used in the fit. For 
the sake of clarity, we rescale the flux via (A/170 /im) . This was 
necessary as there is typically an order of magnitude difference 
between the the far-IR and sub-mm fluxes, and in this way the 
sub-mm, which is the focus of this work, is emphasized. The solid 
line corresponds to /3=1.3, the dotted line to /3=1.5, the short- 
dashed line to (3=1.7, and the long-dashed line to (3=1.9. The 
truncated temperatures shown are in units of -p-j-jy and corre- 
spond to the /3=1.5 and 1.7 cases respectively. See text for an 
explanation of this choice as well as further discussion on the fits. 



table are our Ssso > 3er detections (excluding Nl-001, and 
Nl-002 which are consistent with being at zero redshift). 
They are consistent with z ~ 0.4 — 0.9 ULIGs, based on the 
average luminosity of this sub-sample being above 1O 12 L . 
The likely lowest luminosity (and redshift) source of this set 
is Nl-048, which may be only a Luminous Infrared Galaxy 
(LIG), usually defined as < 1O 12 L.0, depending on the exact 
dust parameters. The bottom half of Table 3 contains five 
LIG candidates (Nl-101, N2-013, Nl-024, Nl-068, and Nl- 
045) based on the range of ((3, Td) we considered. When com- 
paring the redshift ranges obtained using the above rather 
approximate method against the known spectroscopic red- 
shifts, we find fairly good agreement between them. This 
means that the assumptions made were probably reasonable. 
However, the two ULIGs with spectroscopic redshifts have 
measured values at the upper and lower end of the estimated 
ranges, perhaps suggesting a slightly different temperature, 
at least in some cases. 

The only other exceptions here were Nl-008 and Nl-039, 
since, with their z spC c = 0.27, they should have been grouped 
with the others above. There is no reason to suspect intrin- 
sically different dust properties in either case. It is worth 



Source 


z 


log(L) [L ] 


D L [Mpc] 


Nl-078 


0.66-1.00 


12.1- 


-12.6 


3700-6100 


Nl-059 


0.64-0.96 


12.2- 


-12.6 


3500-5900 


Nl-064 


0.58-0.89(0.91) 


12.1- 


-12.6 


3200-5300 


Nl-040 


0.49-0.76(0.45) 


12.0- 


-12.5 


2600-4400 


Nl-048 


0.33-0.56 


11.7- 


-12.3 


1600-3000 


Nl-101 


0.18-0.35 


11.0- 


-11.8 


810-1700 


N2-013 


0.15-0.30 


11.1- 


-11.9 


684-1500 


Nl-024 


0.14-0.30 


11.1- 


-12.0 


621-1400 


Nl-068 


0.08-0.22(0.22) 


10.4- 


-11.5 


330-1000 


Nl-045 


0.08-0.22(0.25) 


10.5- 


-11.6 


320-1000 



a The ranges correspond to different Tj) combinations, the 
first being (1.7,30K) the second (1.5, 40K). The redshifts in 
brackets are measured spectroscopically. The top half of the ta- 
ble has the > 3cr detections, while the bottom half of the table 
has the other sources from the sample with redshifts from the fits 
having lower limits greater than zero. The sources are arranged 
by decreasing redshift. 

noting that there are two radio sources within the Nl-008 
ISO beam. The missing sub-mm flux of Nl-039 may be re- 
lated to the fact that it is one of the sources which were 
solely observed during the poorer weather conditions of the 
May, 2001 observing run - rms = 2.3 mjy (see Table 1), 
compared to a mean rms of 1.5 mjy for the sample. 
Comparing with the redshift estimates derived from the far- 
IR/radio correlations (Section 3.3), we find broad agrement. 
The DCE relation generality gives good overlap with the 
other methods of estimating redshift, while the CY relation 
gives redshifts which are somewhat high. This is because the 
effective temperature in the DCE relation is more appropri- 
ate for our sample. 



3.7 Star Formation Rates 

The first star formation rate (SFR) estimator we consider is 
based on the FIR luminosities. This uses the fact that the 
thermal dust emission is reprocessed stellar light which was 
absorbed primarily in the UV (i.e. where young stars are 
the main contributors). For a discussion of the uncertain- 
ties associated with this and other SFR estimators see e.g. 
Schaerer et al. (1999). Essentially, the problem is to find a 
universal calibration given the vastly varying conditions (e.g. 
different dust properties, and contribution of cirrus, stars, or 
AGN to the thermal spectrum) from galaxy to galaxy. Here 
we adopt a form, which is based on stellar evolution models 
using a Salpeter IMF (Chariot et al. 2002): 

SFRpir 1 n lr> -io^FiR 



[M Q yr- 



= 1.7 x 10" 



L 



(2) 



The second method is based on the radio continuum, where 
the idea is essentially the same as above, but use is also made 
of the well known far-IR/radio correlation (see Section 3.3). 
Here we use an empirical relations based on IRAS galaxies 
with radio observations (Yun et al. 2001): 



SFR rat J 



(5.9 ± 1.8) x 10 



-22 il.4GHz 



[Meyr- 1 ] ™ ' [ WHz 

The radio luminosity is then obtained via 



(3) 
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Table 4. 


Estimating the Star Formation Rates" . 


Source 


z 


SFR^nfM^yr" 1 ! 


SFR FT R [M m yr 


Nl-078 


0.66-1.00 


230-640 


240-730 


N 1-059 


0.64-0.96 


530-1460 


250-740 


Nl-064 


0.58-0.89 


170-460 


200-620 


N 1-040 


0.49-0.76 


160-450 


170-580 


N 1-048 


0.33-0.56 


70-240 


80-310 


Nl-101 


0.18-0.35 


20-80 


20-100 


N2-013 


0.15-0.30 


10-50 


20-150 


Nl-024 


0.14-0.30 


20-110 


20-160 


Nl-068 


0.08-0.22 


3-30 


4-60 


N 1-045 


0.08-0.22 


3-30 


5-70 



a This table is arranged in an same way as Table 4. SFR ra( j uses 
the relation of equation (3), while SFRpjR uses equation (2). 

logLl.4GHz = 20.08 + 21og-D L +logSl.4GHz, (4) 

where L is in WHz -1 , Dl is in Mpc and Si.4GHz is in Jy. 

We present the results in Table 4. The two estimators 
agree with each other reasonably well (as expected for star- 
forming galaxies). The strongest deviation is observed for 
Nl-059. This source has a greater radio flux by a factor of 
roughly ~ 2 from its peers, which may be an indication 
of AGN contribution. In general, estimates for lower lumi- 
nosity sources are more affected by the various systematic 
uncertainties discussed above. The conclusions are: that the 
higher-z candidates have L ~ 10 12 L Q and SFRs of typically 
a few hundred M0yr _1 ; and that the other galaxies with 
estimated redshifts above have L ~ 10 11 Lq and SFRs 
of typically a few tens of M Q yr _1 . The rest of our sample, 
which are more nearby, have lower Lpir and SFRs, consis- 
tent with being more normal star-forming galaxies. 

4 DISCUSSION 

4.1 Summary of spectral properties 

The multiwavelength photometric analysis of the sample of 
galaxies presented in the previous sections provides us with 
an insight into the brightest contributors to the CIB. Our 
sample also holds information on galaxy evolution roughly 
in the range z ~ — 1. 

The series of scatter plots in Fig. 3 showed that a group 
of 5 sources (Nl-040, Nl-048, Nl-059, Nl-064, and Nl- 
078) stand out from the rest in the far-IR/sub-mm, and 
sub-mm/radio projections. Their position in the far-IR/sub- 
mm plot can be explained either by their being colder, or 
at somewhat higher redshift than the rest of the sample 
(the Td/(1 + z) degeneracy). However, when the sub-mm 
slope (S450/S&50) alone is examined, these sources do not 
stand out as one might have expected if their intrinsic SED 
shapes were substantially different from the rest of the sam- 
ple. Thus, we assume an approximately constant SED shape 
across the sample, and arrive at a combined best-fit single 
grey-body (from the S170/S450 and S^so/Ssso slopes) with 
13 ~ 1.5 and T d /(1 + z) ~ 30 K. 

When the sub-mm/radio relation is examined, the same 
group of 5 stands out again, with two redshift estimators 
agreeing with their being at higher (~0. 4-1.0) redshifts. Due 
to uncertainties in the spectral indices, and more impor- 
tantly scatter in intrinsic galaxy properties these indicators 
have Az ~ 0.5. The DCE (Dunne, Clements & Eales 2000) 



relation seems to provide the best-fit to the available spec- 
troscopic redshifts, and may be more reliable. The above 
five sources are the only ones for which this relation gives a 
redshift > 0.4^ 

A different way to look at the data is the sub-mm/near-IR 
relation (Fig. 4) . Here, the segregation of the higher redshift 
candidates is most pronounced - more than 2a separate each 
of the above 5 sources from any of the rest of the sample (al- 
though the location of Nl-032 and Nl-034, being faint at K 
and faint with SCUBA, is poorly understood at this point). 
This projection has the advantage of sampling two spec- 
tral regions with completely different emission mechanisms: 
thermal dust emission vs. stellar light plus dust attenuation. 
This means that the Ta/(1 + z) degeneracy is partially bro- 
ken. This relation would suffer from a completely different 
set of systematic uncertainties than the radio/sub-mm pho- 
tometric redshifts. Such an approach was already discussed 
by Dannerbauer et al. (2002) in the context of their mm- 
selected galaxies compared with the IRAS-selected galaxies 
of the SLUGS sample (Dunne & Eales 2001). However, since 
no robust K-Sgso relation is known at present, we can do 
little beyond obtaining a qualitative confirmation of the rel- 
ative redshifts of our sources. 

Finally, we used the knowledge of the general trends in 
our sample, inferred from the above steps, to attempt to 
constrain some of the properties of the individual galaxies. 
We fitted single, optically-thin, grey-bodies to the 170 /im, 
450 /im, and 850 /im points. For the sample as a whole only 
f3 ~1.5-1.7 provides a good fit, while 1.3 is a poor fit 
for many of the probably local sources, and 1.9 is a poor 
fit for our higher-z candidates. Since the fits only provide 
Td/(l + z), a dust temperature needs to be assumed in order 
to obtain a redshift. We esimate that an acceptable range is 
the (/?, T) combinations (1.7, 30 K)->(1.5, 40 K). These give 
redshifts which are in reasonable agreement with all rela- 
tions examined so far (including the DCE, and CY redshift 
indicators). We used this range to estimate the luminosities 
and SFRs of the 5 high-z candidates, along with some in- 
termediate sources which are possible LIGs up to 2 ~ 0.3, 
representing the tail of the more local bulk of the sample. 
Our results, from near-IR to radio, are consistent with hav- 
ing a sample of mostly local galaxies, some slightly higher 
redshift LIGs and a handful of probable ULIGs at redshifts 
in the range 0.4 < z < 1.0. 

4.2 Bimodality 

Some of the scatter plots discussed above suggest a bi- 
modality in our sample. Whether our particular observa- 
tional selection effects result in bimodality in redshift, or 
just a higher-z tail may be an important point for distin- 
guishing various galaxy evolution models (discussed in the 
next sub-section). Due to the small size of our high-z candi- 
date sample, it is difficult to test their distribution properties 
in detail. However, a simple test which we can perform in- 
volves comparing the \ 2 resulting from fitting a single line 

t An additional candidate is N2-013, but on the basis of the SED 
fits it appears to be a luminous infrared galaxy at perhaps up to 
z ~ 0.3, which is consistent with the redshifts provided by the 
sub-mm/radio relations. 
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Figure 6. Here we test the hypothesis of our sample being bi- 
modal by comparing the \ 2 °f a single-line fit for the entire sample 
(right panel) to a two-line fit for each sub-sample (left panel). The 
dashed lines are ±lcr where a is the rms scatter in the y-dircction. 
Notice that, apart from Nl-048, even with the single-line fit to 
the entire sample, our high-z candidates are > 2a away from the 
best-fit line. See Fig. 3 for representative errorbars. 

(y — mi), against that for fitting two lines, both with zero 
y-intercept. We choose to focus on the £170 vs. Ssso projec- 
tion as here the bimodality is implied, but is not as clean 
as in the sub-mm/near-IR relation. We show the result in 
Fig. 6. We performed this test simply with the ID \ and 
assuming constant error for each source (since they are fairly 
uniformly distributed in any case). The single line fit results 
in x 2 of 108, while the two-line fit results in \ 2 of 37, with 
N—31 (minus either one or two constraints), and for the two- 
line fit we use the nearest line for each point. The two-line 
model thus provides an adequate fit to the data, suggesting 
that each source is drawn from either one population or the 
other. This supports the idea that a handful of our sources 
are a distinct population and lie at z ~ 0.4 — 0.9, while most 
of the sample are at z ~ 0. 

4.3 Comparison with evolutionary models 

The sources studied here are a represenetative sample of 
the brightest ~10% of sources contributing to the CIB. 
They thus provide a test of the various evolutionary mod- 
els abounding in the literature. Models which are consis- 
tent with both the observed CIB intensities, and the number 
counts obtained by various surveys, imply that the majority 
(~80%) of the CIB near its peak (~200 /im) will be resolved 
by sources in the range < z < 1.5. The same redshift range 
sources contribute only ~ 30% of the 850 /im background 
(Elbaz et al. 2002; Chary & Elbaz 2001). Such models re- 
sult in a peak of the SFR density at z ~ 1.0, and then have 
SFR essentially flat until z ~ 4. In general > 70% of the 
star formation takes place in galaxies with Lfir > 10 11 L Q 
(Chary & Elbaz 2001). 



From the redshifts we infer for our sample, it seems to span 
the crucial epoch over which the strongest evolution of the 
SFR density takes place. In general there is no way to fit 
the FIRBACK counts without strong far-IR evolution over 
at least this redshift range. Adopting the starburst tem- 
plate from Lagache et al. (2002), we see that sources less 
luminous than about 10 12 L Q fall below the FIRBACK de- 
tectability beyond redshift ~0.4. Since the FIRBACK se- 
lection excludes normal galaxies beyond z ~ 0.1 and LIGs 
beyond z ~ 0.3, but allows for higher-luminosity sources 
up to z ~ 1.0, our mix of normal, starforming galaxies, in- 
cluding a few possible LIGs, and a handful of most likely 
higher-z ULIGs is in good qualitative agreement with this 
model. The bimodality which this hints at for our selection 
(and which we appear to observe) , is more directly shown by 
a number of specific models (Lagache et al. 2002; Chapman 
et al. 2002; Chary & Elbaz 2001; Dole et al. 1999; Wang & 
Biermann 2000; Franceschini et al. 2001). An easy way to 
achieve such a bi-modal distribution is to phenomenologi- 
cally decompose the luminosity function into a component 
of normal, quiescent galaxies, and a much more luminous 
component of ULIGs, and then have the luminous compo- 
nent evolve more strongly than the quiescent one so that it 
dominates the luminosity funciton by z ~ 1 (Dole et al. 1999; 
Wang & Biermann 2000; Franceschini et al. 2001). This ap- 
proach was exploited in the context of the FIRBACK sources 
(Lagache et al. 2002), showing that a double-peaked redshift 
distribution is predicted, as appears to be confirmed by our 
data. 

Another approach to modelling the N(z) distribution is that 
of Chapman et al. (2002) where the entire infrared luminos- 
ity function is evolved. This combines the colour (i.e. tem- 
perature) distribution of local galaxies with a strong lumi- 
nosity evolution, such as in Xu (2000), to produce the N(z) 
distribution for the FIRBACK population. A bimodallity 
can only be produced here if a bivariate T d distribution is 
used (i.e. by including cold luminous sources), but this is 
strongly evolution dependent. Discriminating in detail be- 
tween such models, including issues such as separating den- 
sity evolution from luminosity evolution, is not currently 
possible. However, once the full redshift distribution of the 
FIRBACK sample is obtained, such discrimination may be 
feasible. 

The FIRBACK selection allows us to investigate the range 
£ z ^ 1. On the other hand z ~ 1 is the lower limit of sub- 
mm/mm selected surveys (Smail et al. 2001; Dannerbauer et 
al. 2002). Thus samples such as ours should act as a bridge 
between the local Universe and the much higher redshift 
population detected in long-wavelength surveys. We illus- 
trated this explicitly in Fig. 4 where we overlaied a number 
of SCUBA-selected sources, and showed that they occupy 
essentially the same sub-mm/near-IR locus as our high-z 
candidates. 

4.4 Conclusions 

We have learned that the brightest ~10% of the CIB is com- 
posed of two different types of galaxy: about 1/6 is the low 
redshift tail of a rapidly evolving ULIG population (sim- 
ilar to the higher-z SCUBA sources); and the other 5/6 
are mainly nearby quiescently star-forming galaxies like the 
Milky Way, with perhaps a few more luminous infra-red 
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galaxies inbetween. This is somewhat contrary to early ex- 
pectations for the nature of the FIRBACK galaxies (e.g. 
Puget et al. 1999) where some sources in the range z ~ 1 — 2 
were expected. This would be the appropriate range for the 
handful of higher- z sources in our sample only if a much 
warmer SED is assumed. However, this is inconsistent with 
the redshift and temperature estimates for the majority of 
our sample, ruling out such models. Further progress on con- 
straining models in detail will come from spectroscopic and 
morphological studies of the entire sample. Understanding 
what makes up the other ~ 90 % of the CIB will need to 
await future far-IR missions with smaller beam-sizes, such 
as BLAST and Herschel, as well as high sensitivity mid-IR 
facilities such as SIRTF. 
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